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\ A number of old and new methods for computing K — * tttt amplitudes on 

Q ' the lattice are reevaluated. They all involve a non-perturbative determination 

. of matching coefficients. I will show how problems related to operator mixing 

£NJ ', can be greatly reduced by introducing the O.P.E. of hadronic currents directly 

on the lattice. Applications to the evaluation of CP-symmetric (AI = 1/2 rule) 
and CP-violating (e'/e) processes are presented. 
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The Kaon sector of the Standard Model may still have some unexpected surprise 
in store for us. Uncovering them requires a very good non-perturbative knowledge of 
hadronic matrix elements. Despite the many non-perturbative approaches developed to 
this aim (chiral lagrangians, 1/iV-expansions,...) PL it seems that only lattice based com- 



p ■ putations can cope with this formidable problem. 

There are two major difficulties, however, which arise in the calculation of hadronic 
matrix elements in lattice QCD. 



1. Decay amplitudes into two or more particles cannot be directly accessed in Euclidean 



space, as a consequence of the Maiani and Testa (MT) no-go theorem ||, except at 
threshold, where, we recall, final state interactions are absent. 

2. Operators can mix with operators of lower dimension with coefficients which diverge 
as inverse powers of the lattice spacing. These contributions must be computed non- 
perturbatively and subtracted, most often leading to prohibitively large statistical 
errors M. 

In this talk I will be mainly concerned with the second of these two problems. In 
Section |I| I recall some of the old methods proposed to deal with the question of operator 
mixing in the case of Wilson fermions (either improved or not) [J In the discussion I will 

* Excellent reviews can be found in refs. [Q-Q- 
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make explicit reference to the problem of explaining the large value of the ratio between 
the 1 = and the / = 2 K — > nn amplitudes (J is the isospin of the two pions in the final 
state). In Section |TTT| I will illustrate a new approach based on the direct use of the O.P.E. 
on the lattice . The idea is to measure in Monte Carlo simulations the small x-behaviour 
(a \x\ <C Aq^D) a lattice spacing) of the hadronic matrix elements of the product of 
two hadronic currents and to compare it with the x-behaviour of its O.P.E. In the O.P.E. 
formula one should look at the Wilson coefficients as known functions (assumed to have 
been computed in perturbation theory), while the matrix elements of the corresponding 
operators are treated as fitting parameters. If the condition a <C \x\ <C Aq^ d is satisfied, 
perturbation theory applies, up to corrections vanishing as some power of a/\x\, and 
the best-fit values of the above parameters will automatically supply the physical matrix 
elements of the finite, renormalized (continuum) operators without the need of knowing 
their expression in terms of bare lattice operators. A feasibility study of this method in 
the case of the two-dimensional 0(3) a-model has been successfully carried out in ref. [ID|. 



Applications of the general philosophy underlying this approach to the lattice computation 
of CP- violating K — > tttt amplitudes, relevant for the evaluation of e'/e fl ]-|J, will be 
presented in Section [TV|. For lack of space I will not discuss here a recent interesting 
extension of the above ideas, in which also the Wilson coefficients are extracted in a non- 
perturbative way from lattice data ||TT|. In this stronger formulation the O.P.E. method 



has been already concretely applied by the authors of ref. |LT| to the computation of the 
first moment of D.I.S. structure functions. Conclusions and an outlook on the perspectives 
of using the lattice O.P.E. method and of its potentialities and limitations can be found 
in Section M. 

II. The old approaches to the AI = 1/2 rule 

One of the major puzzles remaining in hadronic physics is the origin of the so- 
called "A/ = 1/2 rule" in non-leptonic kaon decays. Decays in which isospin changes 
by AI = 1/2 are greatly enhanced over those with AI = 3/2. One finds experimentally 
(neglecting CP violation) 

A(K^nn[I = 0}) ^ 

A(K -> 7T7r[J = 2]) ~ ' y } 

Although the origin of this large enhancement is not theoretically well understood, we do 
know that, in a QCD-based explanation, most of the enhancement must come from long 
distance, non-perturbative physics, or else new physics is winking at us here! 

t For the use of staggered fermions see ref. ||. For the new proposals based on Ginsparg- Wilson, overlap 
or domain wall fermions see the many interesting contributions to these Proceedings. 
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Let me briefly discuss the source of the difficulties in calculating A(K — > tttt) in 
lattice QCD 0. 

For scales below M w , but above the charm quark mass, the AS* = 1 part of the 
effective weak Hamiltonian, T~i£g =1 , can be written in the form 



n 



AS=1 
eff 

<±) 



U V2 

(«Tu^)(«7S)±(s7»(«T^) 



(2) 
(3) 



where 7^ = 7 M (1 — 7s)/2, A u = V^V^,, Gp is the Fermi constant and fx is the subtraction 
point. Throughout this paper I will denote finite renormalized operators with a^ on top. 

The operators and O* - ) have different transformation properties under isospin. 
In particular, is pure I = 1/2, whereas contains parts having both 1 = 1/2 and 
1 = 3/2. An explanation of the AI = 1/2 rule thus requires that the K — > tttt matrix 
element of G_0^ be substantially enhanced compared to that of C + 0^ + \ 

Part of the enhancement is provided by the ratio of Wilson coefficients, G_/G + , 
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their renormalization group evolution from fi ~ M w down to fi ~ 2 GeV, where one finds 
|C_/C + | ~ 2. This factor is, however, too small by an order of magnitude to explain the 
AI = 1/2 rule. The remainder of the enhancement must come from the matrix elements 
of the operators, and these are the quantities that we wish to evaluate on the lattice. 
Attempts in this direction date back to the works of ref. 
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There is a long list of methods proposed in the literature to deal with the problem of 
computing on the lattice the hadronic matrix elements that enter in the non-leptonic weak 
decay amplitude. They are all aimed at bypassing the two major difficulties mentioned in 
the Introduction. In this section I wish to briefly highlight the merits and the drawbacks 
of the most promising among them. Before doing that let me start by recalling what is 
the mixing pattern of the operators O^. 

II- 1. The mixing of on the lattice 

As is well known, the explicit breaking of chiral symmetry due to the presence 
of the Wilson term in the lattice action induces the mixing of the bare operators O^' 
with operators belonging to different chiral representations. Taking into account the 
symmetries left unbroken by the lattice regularization, it can be shown that, in terms 
of bare operators, the renormalized, finite lattice operators possessing the correct chiral 



transformation properties have the form [13 



o(±) = djg + 6 



}(±) 

'PV 



(4) 
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with 

A(±) _ 7 (±)r (±) , V 4 r (±) n (±) , r (±)^ F d , /-(±)^1 
— Z/ PC l^pc z^i=i ^iPC u iPC "t O o-p bv^vF -f- o Id £>aj 

(5) 

OS? = 4?[0$ + C^sa^d + C^ d s l5 d] . 

0^ and Opy are respectively the parity- conserving (PC) and the parity-violating (PV) 
parts of the operators O^l The Oj-p^'s (i=l,...,4) are (6 dimensional) four-quark op- 
erators whose explicit expression can be found in |TJJ]. As we have indicated above, 
parity-conserving and parity-violating operators renormalize separately, because strong 
interactions conserve parity. 

The mixing pattern shown in eqs. @ is the same both with the standard Wilson 
action and with the "tree- level" clover- improved action [ I5| , provided in the latter case 
quark fields are appropriately improved ("rotated"), according to q — > (1 + am /2)q, 
q — > (1 + am /2)q, with m the bare quark mass parameter entering in the fermion 
action. If the action is improved non-perturbatively, i.e. beyond tree- level [|1^], the form 
of eqs. (H) will become much more complicated because of the presence of many more 
operators, needed to correct all higher orders in a, beyond the leading ag^ n log™ a terms 
already taken care of by the tree-level Pauli term in the clover action. To my knowledge 
a complete analysis of the mixing in this case is still lacking. 

Let us now examine in turn the various terms in eqs. fl5|). 

• Dimension 6 operators 

Spin and color structure of the operators of dimension 6 contributing here is the same 



as that of the AS = 2 case discussed in ref. |14|| . Only the flavor structure is different. 
Notice that there is no mixing with dimension 6 operators in the parity- violating part. 
• Dimension 5 operators 

Thanks to the GIM cancellation mechanism, the coefficients Cj?j? and are 

' crF 

actually finite, because the potential 1/a divergence in the matrix elements of and 
Opy gets replaced by a m c — m u factor. Furthermore the CPS symmetry [I7J (CPS = 



CP x symmetry under s — > d exchange) makes proportional to m s — and thus 
vanishing in the limit of exact vector flavor symmetry (exact SU(Nf)y). 
• Dimension 3 operators 

The GIM mechanism softens the a priori possible 1/a 3 divergence of and Cf~? d , 
because, as I said, it has the effect of substituting one factor 1/a with m c — m u . In the 
case of Cgirf an extra 1/a power is taken care of by CPS symmetry which has the effect 
of replacing it by a m s — rrid factor. 



4 



II-2. CP-conserving K — > txtx amplitudes 

In this section I wish to briefly describe a number of methods that have been pro- 
posed to deal with the problem of computing CP-conserving K — > tttt amplitudes on the 
lattice. 



1. In ref. JHJ it was suggested to work with m s = m d and calculate the amplitudes 

(6) 



= (^(^ = 0)^(^ = 0)10^(^)1^ = 0)) 



m s =m d 



with all three particles at rest. We have seen that, setting m s = ntd, GIM cancellation 
and CPS-symmetry cause all mixings in the parity-violating part to vanish, removing the 
need for subtractions altogether. Since one is working at threshold with the two pions 
at rest, the MT no-go theorem || does not apply. So both difficulties mentioned in the 
Introduction are overcome with this simple choice. The method requires, however, a large 
extrapolation from the unphysical, off-shell point, m s = m<j, to the physical one with 
the use of chiral perturbation theory This procedure can be particularly delicate if 



lattice numbers are extracted from quenched simulations, due to the presence of chiral 
logarithms [p0[| . 

2. An alternative method consists in working with the non-perturbatively 0(a) 
improved clover action (for which there are no errors of 0(a) in the spectrum ]I3 and on- 



shell matrix elements of improved currents obey the continuum chiral Ward identities up 
to 0(a 2 ) |T§]). Choosing quark masses such that m K = 2m n , one measures the K — > nn 



amplitudes again with all particles at rest. Since in these kinematical conditions, the 
total momentum transfer vanishes (Ap = 0), one can prove that the matrix element of 
the dangerous, power divergent, subtraction of the operator s^^d is now of 0(a) rather 
than 0(1) and vanishes in the continuum limit. Furthermore pions are at rest and again 
the MT no-go theorem does not apply. Besides the need of a (perhaps less severe) chiral 
extrapolation (unlike the previous case here we are dealing with on-shell amplitudes), 
a problem with this method is the difficulty of tuning quark masses with a sufficiently 
high accuracy to have the condition mx = 2m n satisfied. This problem can be, however, 
alleviated by explicitly performing the subtraction of the s'jsd operator. 

3. The old proposal of ref. [0 is based on the evaluation of the K — > it matrix 



elements of (the positive parity part of) the weak Hamiltonian and makes use of chiral 
perturbation theory, in the form of Soft Pion Theorems (SPT's), to connect K — > 7r to 
K —>■ 7T7T amplitudes. Since only single-particle states are involved, there are no problems 
with the MT no-go theorem. The disadvantage of the method is that, as discussed before, 
the operator mixing problem for the positive parity part of O^ is much more severe than 
for their negative parity part. This makes an accurate evaluation of the matrix elements 
of the renormalized operators very difficult. 
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The relation between K — > txtx and K — > tx amplitudes can be found from the 
classical works of refs. |21| . At leading order in chiral perturbation theory the physical 
amplitude takes the form (for Ap = 0) 

(n+n-\0^^)\K°) = i 7 <±) m *- < . (7) 

Jit 

As the coefficients 7 ( - ± - 1 appear also in the expression for the K — > tx matrix element 

(7T + (p)\d^)\K+(q)) = -6& + 7 (±)p • g , (8) 

by studying this matrix element on the lattice as a function of p ■ q, one can, in principle, 
determine 7^, from which one then obtains the K — > txtx amplitudes 

In order to construct the finite renormalized lattice operators it was suggested 
in reference p2[ to use perturbation theory to determine the finite mixing coefficients 
(CiPC an d C^f) an d subtract non-perturbatively the operator sd, as its mixing coefficient, 

(see the first of eqs. (0)), is quadratically divergent. This approach has been tried in 
ref. 0] with no success as numerical and statistical errors coming from the power divergent 
subtraction completely obscure the physical signal. 



III. Lattice O.P.E. 



In this section I want to illustrate a strategy which, in principle, avoids all the dif- 
ficulties caused by mixing with lower dimension operators and automatically gives the 
physical matrix elements of the effective weak Hamiltonian with the correct normaliza- 
tion 0. 

The method is based on the idea of studying at short distances (a \x\ <^ Aq^ d ) 
the x-behaviour of the O.P.E. of two hadronic currents on the lattice. It does not use chiral 
perturbation theory, and thus in principle applies equally well to the AS = 1, AC = 1 
and AB = 1 parts of the weak Hamiltonian. In addition, it allows one to construct an 
improved weak Hamiltonian (i.e. one having errors of 0(a 2 )), if the improved version of 
the weak hadronic currents |16| is used. The approach is speculative in the sense that it 
is likely to require more computational power than is presently available, although it may 
become practical with the advent of Teraflop machines. 

* For completeness I recall that, in the same notations used in eqs. (0) and (0), one has (0\OW\K°) = 
i{m 2 K - mDS^. 
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I recall that the standard construction of the non-leptonic weak Hamiltonian begins 
with the formula 

n^ = 9 -fj d*xD w (x;M w )T(j pL (x)Jl L (0)) , (9) 

where 

is the longitudinal part of the P^-boson propagator and J p l is the (left-handed) hadronic 
weak current. One then introduces the Wilson operator product expansion in the r.h.s 
of eq. (|9]), a step which is justified by the observation that the dominant contribution to 
the integral comes from very small distances, |x| <C For physical amplitudes, one 

obtains in this way 

(h\nZ\h') = ^YlCi^,M w )M^(h\d^^)\h') , (ii) 

where \h) \h!) are hadronic states, di is the dimension of the operator 0^ l \fx) and Gf/\/2 = 
g 1 w j%M^ /v . The functions Ci(fi,M w ) result from the integration of the Wilson expansion 
coefficients, Ci(x; fi) (defined in eq. (|13D below), with the W^-propagator 

Q(fi,M w )M e w - d > = J d 4 xD w (x;M w ) Ci (x;n) . (12) 

The 0' 8 '(//)'s are quark and/or gluon operators renormalized at the subtraction point /i. 
The functions Ci(fi,M w ) are evaluated in perturbation theory and their running with /i 
is dictated by the renormalization group equation which follows from the /i-independence 
of the l.h.s. of eq. (|TT1) . The sum in the expansion (|TT1) is over operators of increasing 
dimension. As the operator dimension of Ti^ is 6, we will have to consider in the following 
only operators with dimensions di < 6, since the contribution from operators with di > 6 
is suppressed by powers of \jM w . 

All the intricacies and complications of operator mixing in the definition of the finite 
and renormalized operators, O^(fi), come about because the integrals in (|9p and ( |1"2"D are 
extended down to the region of extremely small \x\. The complicated mixing pattern of 
the 0' l '(/i)'s in terms of bare operators arises from contact terms when the separation of 
the two currents goes to zero, i.e. when \x\ is of the order of a. This observation suggests 
that a simple way to circumvent these difficulties is to directly determine the matrix 
elements of renormalized operators by enforcing the validity of the O.P.E. on the lattice 
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for distances \x\ much larger than the lattice spacing a, but much smaller than Aqq D (and 
inverse quark masses), i.e. in a region where perturbation theory determines the form of 
the Wilson expansion. 

We should imagine proceeding in the following way. If J p l is the appropriately 
renormalized (and possibly improved) finite lattice current operator, one starts by mea- 
suring in a Monte Carlo simulation the hadronic matrix element (h\T(J P L(x)Jp L (0))\h'), 
as a function of x in the region a <C \x\ <C Aq^d- The numbers {h\0^> (fi)\h') (eq. flTT|)) 
are extracted by fitting the x-behaviour of (h\T(J P L(x) J pL (0))\h') to the O.P.E. formula 

(/ i |r(j pL (x)J p t L (0))|/ i ') = Ec 4 (x;/i)(/ i |O«( /U )|/ i '), (13) 

i 

where the Wilson coefficients Cj(x;/i) are determined by continuum perturbation theory 
using any renormalization scheme we like. The scale \x should be chosen so that also the 
inequalities a <^ 1/fM <C Aqq D are obeyed. Since we only consider operators of dimension 
6 or lower, the lattice T-product differs from the right-hand side of eq. (|T^) by terms of 
0(|x| 2 Aq CD ), which is then an estimate of the size of the systematic errors intrinsic in this 
procedure. As a last step we insert the fitted numbers {h\0^' (fj)\h') in (|TT|), obtaining in 
this way an explicit expression for the matrix elements of n^- 

The procedure illustrated above requires the existence of a window, a <C <C 
Aq^ d , in which the distance between the two currents is sufficiently small that perturba- 
tion theory can be trusted, but large enough that lattice artifacts, which are suppressed by 
powers of a/\x\, are tiny. For such a window to exist we need to have an adequately small 
lattice spacing. At the same time the physical volume of the lattice must be sufficiently 
large to allow the formation of hadrons. 

A few remarks may be useful at this point. 

• The method determines directly the "physical" matrix elements of the operators 
appearing in the O.P.E. of the two currents, i.e. the matrix elements of the finite, renor- 
malized operators O^'(fi), without any reference to the magnitude of the PU-mass. This 
means that it will not be necessary to probe distances of 0(1/ M w ) with lattice calcula- 
tions. 

• Since it is the continuum O.P.E. which determines the operators appearing in the 



lattice expansion ( ]TT|) , these are restricted by the continuum symmetries. For Aq CD ^> 
\x\ ^> a, the lattice O.P.E. matches, in fact, that of the continuum with discretization 
errors suppressed by powers of a/\x\. 

• Unlike the methods discussed before, this approach automatically yields hadronic 
amplitudes that are properly normalized (in the renormalization scheme in which the 
Wilson coefficients appearing in eq. flT3"|) are computed). 
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As for the applicability of this strategy to the actual case of the CP-conserving AS = 
1 processes, fortunately in the expansion (|Tl"D there appear no operators of dimension 
lower than 6. If lower dimension operators were present they would dominate at short 
distances, since their Wilson coefficients would diverge as powers of 1/x (up to logarithmic 
corrections). In this situation it would be virtually impossible to pick out the matrix 
elements of the interesting dimension 6 operators. 

Operators of dimension 6 have Wilson coefficients which vary logarithmically with 
\x\. At leading order their expression is of the form 

Cl (x; p.) cc [^jf-j ^ = 1 + ^7?° log{\x\n) + . . . , (14) 

where 7^ is the one- loop anomalous dimension of the operator and (3q is the coefficient 
of the one-loop term in the /3-function. 

In the case of the AS = 1 part of Ti.^, the operators which can appear in (|Tl"| ) are 

(eq. (§)) and in addition 

0' = {ml-ml)s{X-%Kd. (15) 

The mass factor in the r.h.s. of eq. (|T5|) comes from a combination of the GIM mechanism, 
which causes O' to vanish when m c = m u , and chiral symmetry, which requires both 
quarks to be left-handed, leading to a GIM factor quadratic in the quark masses. Since 
O' has dimension 6, its coefficient function depends only logarithmically on \x\. 
The anomalous dimensions of the three relevant operators (O^' and O') are 

7^=4, 7i _) = -8, 7o = 16- (16) 

Actually, the contribution of O' to the r.h.s. of eq. flT3|) can be determined separately 
(since its matrix elements do not require any subtraction) and needs not be fitted. The 
anomalous dimensions of the operators O^' are well separated from one another, so it 
might be possible to determine the amplitudes (h\0^> (fx)\h') and obtain the physical 
matrix elements of Ti^a =1 . 



IV. CP-violating AS = 1 processes 



For CP- violating processes in kaon decays (or for B decays), where top-penguin 
diagrams enter at a Cabibbo-allowed level, the strategies described at the points 1. and 
2. of sec. |11-2| for the negative parity operators fail, because the GIM mechanism is not 
operative, as the top quark is too heavy to let it propagate on the lattice. In particular 
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0(±) 

mix with penguin operators (see below). This makes the calculation of the mixing 
matrix of comparable difficulty to that for the positive parity operators discussed at point 
3. of sec. [11- 2j . For positive parity operators, the analysis carried out in sec. [11-1| applies 
with the difference that the mixing coefficients of the color magnetic operator and scalar 
density become more divergent p2j and this will make the numerical determination of the 
renormalized operators less precise. In order to circumvent these problems we propose 
two methods involving a fictitious top quark (with mass m t ) which is momentarily taken 
to be light enough to propagate on the lattice. 



IV- 1. The Renormalization Group method 

The basic idea of what we may call the "Renormalization Group method" is to work 
with two different scales: the first, /i, is larger than m t , so that the corresponding operator 
basis is as in the previous sections; the second, //, is smaller than m t so that a full set of 
new (penguin) operators is generated. The matrix elements of the operators renormalized 
at the scale [i are computed numerically following one of the strategies explained in sec. [n[ 
By matching the result to the amplitude expressed in terms of operators renormalized at 
fi', we extract their matrix elements. In this way, at least in principle, we can obtain the 
matrix elements of the penguin operators without directly computing them. Let me now 
present the details of this procedure. 

At scales a -1 3> M w ^> \i m t , where the fake top quark is active, the CP- 
violating part of the AS" = 1 part of the effective weak hamiltonian, T~t£s =1 \cp, takes the 
form (see eq. (Q)) 

K a / =1 Icp = K^[Ci{toM w ){d\{pL) -6 c M) + 

+C 2 ( f x,M w )(dl( f i)-d c 2 ( f i))}, (17) 



Ol = {srfd)(qrfq) q = u,c,t (18) 
O q 2 = (srfq)(qrfd) q = u,c,t (19) 



where X t = V t dV t * s . The operators O^ 1 ' of eq @ and the corresponding Wilson coefficients, 
C±, in eq (|2|) can be immediately written as linear combinations of the operators and 
coefficients above. 

At scales // below mt, when GIM is not operative, the form of the effective weak 
Hamiltonian becomes 

s-i N 

nt s s=1 \cp = -^J2c j (v',M w ,m t )d j (»') (20) 



V2 



3=1 



10 



where now many new operators appear with a complicated mixing pattern with lower 
dimensional operators |. In fact in the absence of GIM cancellation, in order to write a 
renormalized lattice version of each of the operators Oj it is required to subtract lower 
dimensional operators with appropriate (power divergent) coefficients, and account for 
the mixing with 6 dimensional operators. This is true for both positive and negative 
parity sectors. The coefficient functions appearing in eq. (|20|) have been calculated up to 



non-leading order in perturbation theory in refs. |24|-p9[. 

In order to avoid the difficulties mentioned above we have to keep the GIM mecha- 
nism operative and this is done by having on the lattice a fictitious top quark with mass 
satisfying 

l/a > m t > m c > A QCD . (21) 

To illustrate how the method works in practice let me restrict the discussion below 
to the simpler case of the negative parity operators. One begins by evaluating on the 
lattice at a renormalization scale satisfying fi ^>m t the two matrix elements 

Mii^m) = {h\d c M - 6\{n)\ti) , 2 = 1,2 (22) 

Since GIM is operative, the analysis of point 2. in sec. [II-2| applies (with m u — > m t ) and we 
can define O^/i) — Oj(fi) in terms of bare lattice operators (by appropriately subtracting 
the operator s^^d, as discussed in ref. ||). 

On the other hand, we can also consider the same matrix elements for a renormal- 
ization scale m t ^> fi' ^> m c . Since in this case the GIM mechanism is not at work, all 
the operators which appear in eq. (|20D can contribute and the amplitude Mi(ji,m t ) will 
take the form 

N 

Mi(^mt) = Y,Z^',Li,mt)(h\d j (ii')\ti) . (23) 

3=1 

The rectangular matrix Z^ can be calculated perturbatively by matching the theory with 
(// ^> m t ) and without (// <C m t ) the fictitious top quark Q The results for Z at non- 
leading order can be reconstructed from the works of refs. p4[— 1[29|| . Notice that in absence 
of QCD corrections one would have G\ = 0, Ci = 1. 



§ A convenient basis of operators when QCD corrections are taken into account can be found in the 
papers quoted in ref. P^| . See also refs. |l| ||. 

" This is exactly the method used to calculate the Wilson coefficients in the theory with the physical top 
quark mass, except that in the physical case one must simultaneously integrate out both W boson and 
top quark. Here we are effectively integrating out the W first and the fictitious top quark afterwards. 
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We now choose a fixed value of \J and let m t vary. The matrix elements of interest, 
(h\Oj(fi')\h') , are obtained by fitting the m t dependence of the r.h.s. of eq. (§3|) to 
Aii(fi,m t ), computed numerically as in eq. (|22|), and using the renormalization matrix 
Zij calculated perturbatively. Since the dependence on m t is only logarithmic, this will 



not be an easy job. The procedure is analogous to our use of the x-dependence in sec. |TI 
to find the renormalized matrix elements of operators appearing in the weak Hamiltonian. 
Having determined the numbers (h\Oj(n')\h'), we can insert them into the expression ( ^0|) 
for the effective Hamiltonian. At this point the constraint M w m ( can be removed since 
the Wilson coefficients of the operators appearing in (|i~7|) can be computed perturbatively 
for arbitrary values of m t , including m t = m t . 

As for the errors involved in this procedure, it should be observed that for an accurate 
determination of the matrix elements of Oi([i) (or of Oj (//')) the condition /i ^> m t ^> 
H* Aqcd is not enough. It must also be required that the typical scale, Ahh>, of masses 
and external momenta appearing in the physical process h' —>■ h be much smaller than 
m t . This is because in the matching procedure terms of 0(A/i/ l //m i ) are neglected. 

IV-2. The O.P.E. method 

An alternative method, in the same spirit of the approach followed in sec. |TTTJ, is 
the following. We can avoid the need for any non-perturbative subtraction by separating 
the two currents in x space and having a fictitious propagating top quark. Thus, as in 
sec. |T|, we directly match the x-behaviour of (/i|T(J pi (x) J^ i (0))|/i')^ 1 p =1 (the subscript 



indicates the presence of the fictitious top) to the formula 

(h\T(J pL (x)JU(0)W)tt 1 = J2 c ^^'^t)(h\d^')\h') , (24) 



j 

where 

J p l(x)4 l (0)\ AS=1 = s(xh L Ax)t(0)l L AO) - s(x)rfc(x)c(0)l L AO) . (25) 

The coefficients Ci j2 (a;; fi',m t ) = Ci j2 (x; fi') are the same as those of sec. |T[ The others are 
complicated functions of the anomalous dimension matrix which can be worked out from 
the results of refs. [^(J and |28j and computed numerically. 

Both methods proposed in this section require small enough lattice spacings to 
accommodate a number of scales. Like the method of sec. |ITJ their full implementation is 
likely to require the next generation of supercomputers. 

V. Conclusions 



The new experimental results on e'/e BO], combined with the old ones [511], have 



prompted a wage of theoretical work aimed at testing the validity of the Standard Model in 
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this crucial corner of the theory. Lattice, as shown recently by the investigation carried out 



in ref . |32] , has the potentiality of providing a clearcut answer to the question on whether 
CP-violation can be understood and described within the framework of the Standard 
Model, as it is formulated today. In my opinion all possible efforts should be addressed 
by the lattice community to this fundamental issue. 

With this in mind in this talk I have reviewed a number of old and new approaches 
aimed at studying the AI = 1/2 rule on the lattice, using Wilson-like fermions (similar 
methods could also be used for staggered fermions ||) and I have discussed a new approach 
which can be equally well applied to CP-conserving and CP-violating AS = 1 processes. 

In particular the methods of refs. |I7fl and p2f involving K — > txtt and K — > 7r 



amplitudes, respectively have been reevaluated. The last approach is likely to be more 
difficult than the first one, because of the large number of mixing coefficients which have 
to be determined non-perturbatively. It may however give complementary information to 
the results obtained with the K — > mi methods (points 1. and 2. of sec. [11- 2| ) and provide 
a check of the accuracy of chiral extrapolations. 

To overcome the difficulties inherent in the construction of finite renormalized lattice 
operators, a new strategy, based on the study of the short distance behaviour of the O.P.E. 
of two hadronic weak currents on the lattice, was proposed in ref. Q. The approach is 
theoretically very appealing and can be applied to both CP-conserving and CP-violating 
AS = 1 processes. There are also good indications |TT| that the approach can be extended 
in a way as to allow a fully non-perturbative evaluation of the (first few) moments of D.I.S. 
structure functions 0. 

An other interesting feasibility study in the direction of testing the O.P.E. method 



in actual simulations has been undertaken in ref. (JTOj] . It consists in the study of the 
small x-behaviour of the O.P.E. in the two-dimensional lattice 0(3) cr-model. A prelimi- 
nary analysis of Monte Carlo data indicates that the measured small x-behaviour of the 
one-particle matrix elements of products of operators matches the logarithmic behaviour 
expected from perturbative calculations, thus allowing the (non-perturbative) evaluation 
of the matrix elements of the operators appearing in the O.P.E. 

The advent of Teraflop Supercomputers may render the O.P.E. method a viable 
strategy also for QCD, allowing us to directly extract physical amplitudes from Monte 



" The generalization proposed in ref. fll|| makes use of (gauge non-invariant) quark and gluon states in 
order to get a sufficiently large set of equations and be able to extract all the unknown coefficients and 
matrix elements. In these circumstances non-BRST operators vanishing by the equations of motion 
can contribute to the O.P.E. |3!j. If they are not included, care must be exerced in choosing the 
external states in order to avoid contaminations from these unwanted operators onto the relevant 
matrix elements. 
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Carlo data, which may compete with the more standard approaches illustrated in sec. [Tl 
or with those, based on Ginsparg- Wilson, overlap or domain wall fermions |34]], that are 
now starting to be developed. 
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